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The Detection of Low Mass Companions in Hyades Cluster 

Spectroscopic Binary Stars 

Chad F. Bender^'^ and Michal Simon^ 
ABSTRACT 

We have observed a large sample of spectroscopic binary stars in the Hyades 
Cluster, using high resolution infrared spectroscopy to detect low mass compan- 
ions. We combine our double-lined infrared measurements with well constrained 
orbital parameters from visible light single-lined observations to derive dynamical 
mass ratios. Using these results, along with photometry and theoretical mass- 
luminosity relationships, we estimate the masses of the individual components in 
our binaries. In this paper we present double-lined solutions for 25 binaries in 
our sample, with mass ratios from ~ 0.1 — 0.8. This corresponds to secondary 
masses as small as ~ 0.15 Mq. We include here our preliminary detection of the 
companion to vB 142, with a very small mass ratio of g = 0.06 ± 0.04; this indi- 
cates that the companion may be a brown dwarf. This paper is an initial step in 
a program to produce distributions of mass ratio and secondary mass for Hyades 
cluster binaries with a wide range of periods, in order to better understand bi- 
nary star formation. As such, our emphasis is on measuring these distributions, 
not on measuring precise orbital parameters for individual binaries. 

Subject headings: binaries: spectroscopic — open clusters and associations: indi- 
vidual (Hyades) — stars: fundamental parameters — techniques: spectroscopic 

1. Introduction 

Observations of spectroscopic binary stars provide dynamical measurements of stellar 
mass and binary mass ratio that are important inputs to theories of binary star formation 
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( iBonnell et al.ll2003l : lTohlindl2002l : IClarkdl200ll ). Measuring the radial velocity versus orbital 
phase for a single-lined spectroscopic binary (SBl) yields orbital parameters and the mass 
function, /(M), in which the mass ratio, q = M2/ Mi, is inseparable from the orbital inclina- 
tion. For an ensem ble of SB Is, a statistical a pproach can be used to derive the distribution 
of mass ratios (e.g.. lMazeh fc Goldberg]|l992l ). If, however, a binary is observed as a double- 
lined system (S B2), then the dyn amical mass ratio follows directly from the radial velocity 
measurements (IMazeh et al.ll2002l ). With an estimate of the primary mass, say based on the 
spectral type, the mass ratio gives the secondary mass, and in the case of an ensemble of 
SB2s, the secondary mass distribution. 

Obtaining a purely dynamical mass ratio distribution using SB2s is challenging: a binary 
with a mass ratio of much less than one has a small flux ratio, a = F2/F1, so the secondary 
component is difficult to detect. This flux ratio problem is particularly inhibitive at visible 
wavelengths where much of the long term monitoring of spectroscopic binaries has occurred. 
For this reason, most identified SBs are SBls, and there is a strong selection effect favoring 
the detection of SB2s with q near 1. A binary composed of main sequence stars with unequal 
masses will have an a that increases towards longer wavelengths, making the companion 
easier to detect using infrared spectroscopy. Su ch observations have allowed the mea surement 
of SB2s with q's as small as ~ 0.1 - 0.2 fe.g.. lMazeh et al.lbood : IPrato et al.ll2002h . 



Duquennoy fc Mayorl (jl99ll ) examined a sample of 164 nearby SBls, SB2s, visual bina- 



ries, and common proper motion pairs, with F- and G-type primary stars, and found that 
the ma ss ratio distribut i on for medium to long period binaries increases towards small mass 



ratios. IGoldberg et al.l (120031 ) used a well defined sample of 129 SBls and SB2s from the 
Carney-Latham sample of high proper motion stars (jCarney et al.lll994l ) to derive a bimodal 
mass ratio distribution with peaks at g ~ 0.2 and q ~ 0.8. They found somewhat different 
distributions for halo versus disk stars, and for bi naries with p r imary masses greater than 
and less than 0.67 Mq. By working in the infrared, iMazeh et al.l (120031 ) observed 32 binaries 
from the Carney-Latham sample as SB2s, with mass ratios as small as ~0.20. When com- 
bined with SB2s from visible light spectroscopy, they found a mass ratio distribution that is 
approximately flat from q ~ 0.3 — 1.0. 



The s a mples studied by iGoldberg et al.l (120031 ) . iMazeh et al.l (120031 ) . and, for g's near 
1 by iLucyl (120061 ). represent averages of star formation outcomes in the solar neighborhood 
over billions of years. Most star formation occurs in localized regions within molecular 
clouds. Studies of the binary population in diverse star forming regions can isolate the 
physical parameters that determine the properties of the binaries formed. For example, 
the competi ng effects of fragment a tion, dynamical interactions, and accretion are not full y 
understood (iGoodwin et al.l 120071 : iBallesteros-Paredes et al.l l2007l : IWhitworth et al.l 120071 ). 
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Observations of binaries in open clusters can clarify the important processes because they 
share a common star formation origin and history. Ideally, such observations would be carried 
out on young clusters that have an easily identifiable membership and still retain their most 
massive members. However, it is difficult to measure the orbital velocities of the youngest 
stars because they tend to have high rotation velocities. Stars in older open clusters still 
retain much of their cluster identity yet are rotationally spun down. 

The Hyades is one such nearby open cluster: it is well studied and presents an excellent 
laboratory for measuring a well defined sample of binary stars in order to inves tigate binary 
star formation. The cluster has an age of ~ 650 Myr (ILebreton et al.l 1200 if), old en ough 
that stars later than spectral type F3 h ave slowed to ysini < 25kms~^ (jKraftlll965l ). and 
G-type and later to vsini < lOkms"^ ( IStauffer et al.l 119871 ). Most mid K-type and earlier 
cluster n aembers have individu al Hipparcos distances, and with a mean cluster distance of 
~ 46 pc (IPerrvman et al.lll998r) the known binaries are bright. Additionally, with metallicity 
[Fe/H]=0.14 (ILebreton et al.l 1200 ll ) the spectral lines of Hyades members are deep and well 
suited to spectroscopic analysis. 

196^ 



Radial velocity surveys of the Hyades (e.g. 



Wilson 



19841 : Istefanik fc Latham! Il985] : ICriffin fc Gum]|l978l . 



194 



1981 



KraftI 



iGriffin et al. 



Detweiler et al. 



19851. 



198 



^ have 



identif ied a sufficient numbe r of SBs to facilitate a statistical analysis of their physical prop- 
erties. iPatience et al.l (119981 ) carried out a study of Hyades "visual binaries" using speckle in- 



terferometry. These wide binaries have periods longer than most SBs and seve ral are the wide 



compo nents of hierarchical triples with known SBs. From their observations, iPatience et al. 



( 119981 ) derived a photometric mass ratio distribution for long period systems that appears 
to increase towards small q. 

At the Harvard-Smithsonian Center for Astrophysics (CfA), Robert P. Stefanik (RPS) 
and David W. Latham (DWL ) have monitored Hyades members for SBs by visi ble light 



spectroscopy since 1979 (e.g.. iLatham &: StefanikI Il982l : IStefanik fc Lathamlll985l ). Their 



observations have refined the cluster membership and yielded precise parameters for many 
SBls and SB2s. In 2003, the current authors began a collaboration with the CfA group to 
extend, through infrared observations, the detection of binary companions to smaller masses 
than possible by visible light spectroscopy alone. To this end, RPS and DWL made available 
to us their unpublished parameters for many Hyades SBls. In this paper we present infrared 
SB2 detections that span the range of mass ratios from g ~ 0.1 to g ~ 0.9. We verified 



these results using the saniple of Hyades SBls available in the literature (e.g.. ISanfordlll921 



pie ( 

Griffin fc Gunnlll978l . ll98ll : lGriffin et al.lll985f ). When combined with available surveys (both 
present and future) of Hyades binaries, our SB2 detections will enable better determinations 
of the distribution of secondary masses. In ^ we present our Hyades binary sample; ^ 
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describes our infrared observations and velocity measurement techniques; §1] includes the 
SBl and SB2 solutions for our sample binaries; ^discusses the quality of our SB2 solutions 
and derives the secondary masses; and §S] presents a brief summary and comments on our 
future papers. 



The Infrared SB2 Sample 



Because the Hyades cluster is nearby, it covers a large angular extent on the sky, so 
determining the cluster membership is difficult. Considerable effort has been applied to 
identify true Hyades members using meas urements of proper motion, radial velocity, parallax 
and color-magnitude relationships (e.g.. 



Griffin et al.lll988l : lReidlll992l : IPerrvman et al 



van Bueren 



1998 



1952 



van Altena 



de Bruijne et al 



1969; Hanson 1975 



20011 ). Nevertheless, 



these surveys contain many candidates that are non-members or have uncertain membership 
status. In particular the membership of SBs is often ambiguous until their center-of-mass 
velocities can be calculated from orbital solutions. The multi-decade long campaign at Cf A to 
monitor candidate Hyades binaries and multiple systems has, for many systems, resulted 
in precise measurements of both the orbital parameters (the period, P, the eccentricity, e, 
the time of periastron passage, Tq, the longitude of periastron, u, the semi-major amplitude 
of the primary, Ki, and the center-of-mass velocity, 7) and the mass function, f{M). 

In 2003, we used the preliminary CfA SBl solutions to select a sample of 32 binaries 
suitable for further observations in the infrared as SB2s. We will refer to this as the infrared 
sample, to distinguish from the entire CfA sample. All infrared sample members are con- 
firmed as cluster members through the methods listed above. To avoid early type stars with 
large rotational velocities (§T]) we restricted the sample to targets with primaries of spectral 
type F and later. Despite this, however, five targets do have Vrot ^ 25 km and required a 
modified analysis ( §3.4p . Most of the infrared observations were obtained with the CSHELL 
spectrometer at the IRTF, which, for Hyades members, has an effective magnitude limit of 
i/~9 in a 90 minute integration; this corresponds to targets with early-M primary stars, and 
serves as the faint cutoff for our sample. We intended the infrared observing campaign to 
last at most a few years, during which the orbital phase of systems with long periods would 
not change significantly. Using f{M) from the SBl solutions, an estimate of the primary 
mass based on spectral type, and letting i = 7r/2, we derived the minimum mass ratio, 
^min, allowed for each binary. From q^i^ we predicted the velocity separation of long period 
systems, and excluded those with separations of only a few kms~^ or smaller. Long period 
systems with large predicted velocity separations were retained, with the expectation that 
observations at multiple phases would not be possible. Table [1] fists the infrared sample. 
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Columns 1-5 give the name and HD number of the primary star, the J2000 coordinates, and 
the primary spectral type as listed in the SIMBAD database and confirmed with the analysis 
described in ^ Columns 6 and 7 give V magnitude {Hipparcos where available and SIM- 
BAD otherwise) and H magnitude (2MASS), respectively. Figure [T] shows the distribution 
of i7-magnitude, from 2MASS, for the infrared sample. The hashed regions indicate sample 
members for which we did not detect a companion, (see §5.2p . 



3. Infrared Spectroscopy and Double-Lined Measurements 

3.1. Observations 

We observed the infrared sample at the NASA Infrar ed Telescope Facilit y (IRTF) with 
CSHELL, the facility near-infrared echelle spectrograph ( Greene et al. 19931 ): Table [2] con- 



tains a log of the observations. CSHELL uses a 256x256 pixel InSB detector with a plate scale 
of 0'.'2 pixel"^, and operates under natural seeing at the IRTF. Echelle orders are isolated 
with an order sorting circular variable filter. All of the observations used a single grat- 
ing setting centered at 1.5548 /xm because this spectral region contains several deep lines, 
a requirement for precisely measuring radial velocities. We used a 0'.'5 slit which provided 
a spectral resolution of ~ 30,000. The spectra we obtained have a free spectral range of 
~ 40A. 

We observed each binary as a series of short integrations, from 60 to 300 s in duration, 
nodding the telescope ~ 10" along the slit between integrations in an ABBA pattern. By 
differencing each set of A and B frames we removed contributions from detector bias, dark 
current, and sky background. We took a series of flat and dark frames at the beginning of each 
observing session, which were median filtered and normalized to provide a fiat field correction 
for each target "A-B" image. To measure the dispersion solution, we observed Ar and Kr 
arc-lamp spectra from CSHELL 's internal lamps. We extracted each "A-B" image using a 
custom optimal extractioii code written in IDL, patterned after the algorithms described by 



Piskunov fc Valentil (120021 ). and averaged the set of resulting spectra for each target. We 
adjusted the total integration time for each target to provide spectra with a signal-to-noise of 
S/N ~ 100, although we relaxed this condition for targets with a large qmin, and also for the 
faintest targets to avoid very long integrations that were an inefficient use of telescope time. 
Observations in October 2003 suffered from poor weather conditions and did not result in any 
usable spectra. During the 2005 and 2006 observations the CSHELL detector experienced 
brief, uncontrolled temperature warm-ups that increased the thermal background by as much 
as two orders of magnitude over the nominal few tens of ADUs, affecting ~ 20% of the 
integrations. Spectra obtained under these circumstances have decreased S/N, but most are 
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still usable. 



The CSHELL observations from 2004 suggested that high S/N observations of H69 and 
vB 142 with a larger telescope would be particularly useful. H69, with if ~ 9.1 mag, is 
a faint target for the IRTF. The SBl solution for vB 142 yields a mass function with a 
very small gmin of 0.04, and the 2004 October 1 SB2 measurement indicated that q was 
indeed small and the companion possibly a brown dwarf. We subsequently observed these 
binaries on 2005 February 2 2 at the W. M. Kec k Observatory with NIRSPEC, the facility 
near-infrared spectrograph (IMcLean et al.ll2000l ). We also observed vB 43, which at that 
time was nearing a phase with a large velocity separation. These observations are included 
in Table [2l We used NIRSPEC in its cross-dispersed echelle mode with the 2-pixel slit, 
which provided a spectral resolution of ~ 31,000. Each target was observed as a sequence of 
ABBA integrations, using a nod of ~ 10". A sequence of flat and dark frames was obtained 
at the beginning of the night and raedian filtered. Night sky OH emission lines, identified 
from the catalog of iRousselot et al.l (120001 ) . provided a simultaneous wavelength reference 
to determine the dispersion solution. We extracted the spectra in IDL using REDSPECEI, 
the NIRSPEC facility software package. In addition to normal image processing, REDSPEC 
rectifies the images in both the spatial and spectral dimensions to extract each order. The 
extracted nod pairs were averaged for each target. We worked only with NIRSPEC order 49, 
centered at ~ 1.56 /xm and spanning ~ 220A, because it is nearly completely free of telluric 
absorption; the other orders pr ovided by our grating setting are contaminated to varying 
degrees (see lBender et al.ll2005l . Fig. 1). 



3.2. Template Spectra 



We measure the radial velocities of a binary spectrum with reference to template stellar 
spectra having spectral type, metallicity, and rotational velocity similar to the binary com- 
ponents. We have two libraries of observed single star tem plates. The first wa s obtained with 
NIRSPEC and covers spectral typ es from GO through M 9 (IBender et al.ll2005l . Fig 3-6). The 
template velocities are tied to the iNidever et al.l (120021) reference frarn e; this procedure and 
the individual template velocities are reported in ISimon et al.l (120061 ). The s econd library 



was obtained with CSHELL and covers spectral types from GO through Ml (IMazeh et al. 



2002 



Fig 2). This library is considerably more sparse in spectral type compared with the 
NIRSPEC library and some CSHELL templates have a low S/N. Both template libraries 
are comprised of stars with small rotational velocities. We used a non-linear limb darkening 



http: / /www2. keck. hawaii.edu/inst /nirspec/redspec/ 



- 7- 



model (IClaretl l2000f) to rotationally broaden the templates as needed, following the proce- 
dure outlined by iGrayl (119921 ) . In the analysis presented here, we relied primarily on the 
NIRSPEC templates because of their high quality, but also utilized the CSHELL templates 
where they matched the spectral types of the targets. 



3.3. Radial Velocity Measurements 
We analyzed the infrared spectra wi th a two-dimensional co rrelation algorithm similar 



to the TODCOR algorithm described by IZucker fc Mazehl (jl994l ). For each infrared obser 



vation, we used the SBl parameters to calculate the orbital phase and the radial velocity 
of the primary. We then estimated the range of possible secondary velocities at that phase, 
corresponding to q ranging from q^i^ to 1, with some allowance for the uncertainty in the 
primary mass estimate (^. We analyzed each target with two or three templates that 
matched the known characteristics of the primary and a range of four to eight plausible 
secondary templates, while leaving the flux ratio unconstrained. This resulted in sets of 
measured velocities and flux ratios for each combination of templates, which we averaged to 
obtain the final measured velocities reported in Table [2j We did not detect the companion 
for seven of the binaries in our sample, and for three others were successful for only part of 
our observations. We indicate these observations in Table [2] as missing data, and address 
these systems in detail in §5.21 

We estimated the uncertainties of both the primary and secondary velocity measure- 
ments as a combination of two factors: (1) the velocity uncertainty of the templates and (2) 
the uncertainty in measuring the center of the correlation peak due to noise in the target 
spectrum. We took the contribution from the templates simply as the variance of the veloc- 
ities measured over the set of template pairs. We measured the contribution from noise in 
each target spectrum as follows. First, we created a model binary from a pair of templates 
used in the target spectrum analysis, with the same velocity separation and flux ratio as 
measured in the target. Then, we added random noise to the model and analyzed it using 
the correlation algorithm with the same pair of templates. We adjusted the amplitude of 
the noise until the peak correlation from the model equaled the peak correlation measured 
in the target spectrum. We then generated 100 copies of this model, each using a unique 
noise vector at the prescribed amplitude. Analyzing these models resulted in distributions 
of primary and secondary velocities, which we fit as normal distributions. We repeated this 
with four different template pairs, and took the average of the distribution widths to be 
the uncertainty contributed by noise in that target spectrum. By combining this average in 
quadrature with the uncertainty due to the templates, we arrived at our best estimate of the 
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uncertainty for an individual velocity measurement. We used this procedure to empirically 
determine the uncertainties on both our primary and secondary velocity measurements for 
each of our observed target spectra. 

The relative contribution to the total uncertainties from (1) and (2) varied for each 
target spectrum. Templates with later spectral type have a worse velocity precision that 



those with an earlier spectral type (IBender II2006I ). So, for example, an analysis using M-star 



templates would have a large portion of the total uncertainty contributed by the templates. 
Conversely, a target spectrum with low S/N, often the case for the fainter binaries in our 
sample, would derive most of its total uncertainty from the measurement of the correlation 
peak. An addition al source of noise that has effected previous correlation analyses (e.g.. 



Bender et al.ll2005l ) comes from a mismatch between the target primary and the template 
primary. If these spectra have different physical properties, particularly the metallicity of 
various species, then the precision of both the measured velocities and flux ratio is effected. 
We examined this on several of our target spectra where the velocity uncertainties obtained 
using the above procedure seemed inadequate (see We reanalyzed the model binaries 
from (2) using alternate template spectra of similar spectral type, thereby simulating po- 
tential mismatch between the templates and the target spectra. We found that the increase 
the velocity uncertainty over that obtained in (2) was negligible. This is likely due to the 
small spectral range of our CSHELL templates, which limits the possibility of metallicity 
mismatch. Consequently, the velocity uncertainties reported in Table [2] do not include any 
contribution from this effect. 

To verify that our visible and infrared velocity reference frames were in good agreement, 
we compared the primary velocities measured in the infrared spectra with those predicted 
from the visible light SBl solutions. Figure [2] shows the distribution of the difference be- 
tween the predicted velocities and the measured velocities. The distribution has a mean 
of ~ 0.3kms~^ and is fit by a Gaussian wi th a width of ~ .9kms~^. This difference is 



smaller than our infrared velocity precision (jSimon et al.ll2006l ) and demonstrates that the 



two reference frames are indeed consistent to within the measurement uncertainties. The 
samples in Figure [2] with a velocity difference larger than ~ 5kms~^ are attributable to 
infrared spectra with a low S/N. 



3.4. Rapid Rotators 

The targets vB 8, vB 30, vB 68, vB 77, and BD+02 1102 have F-type primary stars 
that are rapidly rotating, with Vrotsinz from ~ 40 — lOOkms"^, and so the spectral lines of 
these stars are strongly rotationally broadened. The limited spectral range of the CSHELL 
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spectra did not contain enough strong features to measure the primary velocities in the 
infrared spectra of these binaries. We considered, however, that lower mass companions 
will have spun down and might therefore be sensitive to a correlation analysis. To examine 
this we modified our correlation analysis to use an artificial primary template comprised of 
a featureless, flat spectrum. We analyzed each target spectrum with a range of plausible 
secondary templates, allowing the flux ratio to vary, to measure the secondary velocity. This 
procedure was successful for vB 68, vB 77, and BD+02 1102, and the velocities we measured 
for these systems are included in Table [21 along with a note indicating the altered analysis. 
The associated uncertainties were modeled using the procedure described in §3.31 and the 
artificial primary template described above. We did not detect the secondaries for vB 8 or 
vB 30; we discuss both systems in §5.2[ 



3.5. Triple Systems 



Four binaries in our infrared sample, L20, vB 102, vB 151, and vB 40, are the inner 
pair of hierarchical triple systems. Speckle imaging from [Patience et al.l (119981 ) identified the 
wide companions in these systems with orbital periods from a few tens to a few hundreds of 
years, assuming circular orbits. We briefly describe how the presence of a third spectrum in 
our observations affected the analysis. 

The wide companion to L20 has similar brightness to the primary of the inner binary. 
Its spectrum is clearly apparent in our infrared observations, with a radial velocity of ~ 
32kms~^. When left unaccounted for, this contribution limited the precision of our radial 
velocity and flux ratio measurements for the inner binary. To correct for this, we used a one- 
dimensional correlation to determine that the wide companion spectrum was well matched 
with a K5 type template and contributed about half of the total flux. We then used the 
template spectrum to subtract out this contribution prior to the normal two-dimensional 
cross- correlation step. In this manner, we increased the precision of our inner binary velocity 
measurements and decreased the uncertainty on K2 by a factor of two over that obtained 
when not accounting for the wide companion. 



Using the component masses derived by [Patience et al.l (119981 ) and the [Baraffe et al. 



( 1l998[ . BCAH) stellar models, we estimate that the wide companion to vB 102 contributes 
only ~ 10% of the total flux at 1.6 /im. Additionally, its position angle ([Patience et al. 



1998[) placed it near the edge or outside of the CSHELL slit during our observations. If the 
outer binary is in a circular orbit, its period is ~ 30 years, and the velocity separation from 
the inner binary 7- velocity is only a few kms~^. Additionally, the radial velocity of the 
inner binary primary was only a few kms~^ from the 7- velocity during our observations. 
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Consequently, any contamination in our spectra from the outer companion was masked by 
the primary and we were unable to detect it; we therefore solved vB 102 as a normal SB2. 

The wide companions of vB 151 and vB 40 have large angular separations and were 
at position angles such that they fell outside of the CSHELL slit during our observations. 
Even had they fallen within the slit, their contributions to the total flux are small at 1.6 /im, 
~ 10% and ~ 2%, respectively. We did not detect any contribution from the wide companion 
in either system, and concluded that our SB2 correlation analyses of the inner binaries were 
unaffected. 



With the permission of RPS and DWL, we list in columns 2-7 of Table |3] the CfA 
orbital parameters for the 25 SBls that we succeeded in turning into SB2s with infrared 
detections; the seven sample members for which we did not detect the secondary ( §3.31) 
are not included. The CfA parameters are based on preliminary orbital solutions that may 
change slightly when the final results are published, but these changes are expected to be 
insignificant for our present purposes because the uncertainties in our final mass ratios are 
dominated by the infrared velocities for the secondaries. 

We have from one to three infrared spectra for most our binaries, and none were observed 
on more than five occasions. This contrasts with the several tens or more visible light 
observations for each system. Consequently, the primary velocities that we measured in 
the infrared do not further improve the precision of the SBl parameters. To solve each 
system as an SB2, we took, without modification, the SBl parameters from Table [31 and 
used a least-squares fitting routine to solve the infrared V2 measurements for K2. For several 
of our binaries (vB 62, vB 69, H532, vB102, vB 142, vB 121, and vB 151), this procedure 
yielded reduced chi-squared, xl^ much greater than one, indicating that the corresponding V2 
uncertainties (§3.3) are underestimated. Assuming that the missing uncertainty contribution 
in V2 is normally distributed, we account for it by scaling the measured K2 uncertainty by 



\rx^. Because we do not understand the source of this additional velocity uncertainty, we 

choose to retain in Table [2] the uncertainties derived in §3.31 To maintain consistent results 
throughout our sample, we applied this correction to all of the binaries; column 8 of Table [3] 
lists the derived K2 along with the adjusted uncertainties. 



4. 



Infrared SB2 Solutions 



4.1. 



CfA SBl Parameters 




We list in column 9 of Table [3] the mass ratio for each system, calculated from Ki and K-^ 
as g = K1/K2. We verified q for the seven SB2s with large xt using the technique of I Wilson 
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( ll94ll ): in each case, these results agree to within la of those from the least-squares fitting. 
The uncertainty we report for q includes the unlisted uncertainty on Ki from the preliminary 
CfA SBl solution. This contribution, however, is generally negligible when compared with 
the large uncertainties on K2. For example, Ki for vB 9 is determined to ~ 6%, which 
contributes only ~ 10% of the total uncertainty on q; most of the infrared sample members 
have Ki determined much more precisely. L20 has K2 determined to better than ~ 4% and 
is one of our most precisely measured infrared systems; Ki here contributes only ~ 1% to the 
total uncertainty on q. Columns 10-13 in Table [3] give the semi-major axes and component 
masses, combined with the unknown sinz. 

Figures [3] and m show the SB2 velocity curves and measured secondary velocities, plotted 
against orbital phase for the binaries in Table [31 The primary velocity curves shown come 
directly from the SBl solutions, while the only free parameter in the secondary velocity 
curve is K2. We show the secondary velocity uncertainties listed in Table |2l without the 
scahng correction described earlier, because these values represent our best understanding of 
the measurements. As a consequence, some of the velocities are shown with underestimated 
uncertainties. Figure [5] shows two different representations of the mass ratio distribution 
for these binaries: (a) shows the case where each q is determined with equal precision; (b) 
distributes each measured q over its corresponding uncertainty. The similarity of the distri- 
butions, within the reported \/N uncertainties, suggests that our underestimate of the V2 
uncertainties has a small effect on the overall distribution. These distributions are obviously 
not complete for the cluster, most importantly because it does not include the large sample 
of binaries with q > 0.6 that are detected as SB2s in the CfA visible light spectroscopy. Fig- 
ure [5] does, however, clearly demonstrate the applicability of infrared observations to binaries 
with small mass-ratios. 



4.2. Alternate SBl Parameters 



Alternate SBl orbits are av ailable in the literature for nine of our infrared bi naries: eight 
by R. F . Griffi n and colleagues (ICriffin fc Gunnlll978l . Il98ll : ICriffin et al.l 119851 ) and one by 
SanfordI (jl92ll ). The parameters published by these authors provide additional verification 
of our assertion that the uncertainties we report for K2 and q are dominated by our infrared 
measurements of V2. For each of these nine systems we used the alternate SBl parameters 
without modification to derive K2 and q, following exactly the procedure describe above for 
the CfA parameters. However, we did not attempt to reconcile th e various velocity zer o 
points, which have small differences of order lkms~^ or less (e.g.. iDetweiler et al.lll984j ). 
Table m lists these SBl parameters, their references, and our alternate SB2 results. 
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The binaries L20, vB 43, vB 62, vB 69, vB 77, H509, and vB 121 have q hsted in 
Tables [3 ] and HI that agree to ~ la or better. The parameters Tq and uj given for vB 40 by 



SanfordI (Il92ll ) have a 180° phase difference from the CfA parameters, which does not affect 
the derivation of q. The shghtly greater than la difference in the derived q for this binary 
could result from differences in the velocity zero point; alternatively, vB 40 is a triple, and 
orbital motion of the wide pair over the ~ 90 years between Sanford's measurement and our 
own could account for the discrepancy. The q calculated for L57 using the Griffin and CfA 
SBl parameters demonstrate the importance of a consistent velocity reference frame when 
combining SBl and SB2 observations ( §3.3p . Due to its long period, L57 has a small Ki, 
which amplifies the effect of any offset in the ve locity reference frames. Solving this system 
using the SBl parameters of iGriffin et al. Jl985h . but with the CfA 7, yields q = 0.85 ±0.24, 



nearly identical to the result in Table [31 This issue is independent of the large uncertainties 
we report on q for L57, which result from our observing it only once in the infrared. 

The results shown in Table H] confirm that the uncertainties in our reported SB2 pa- 
rameters originate primarily from our infrared measurements, not from the underlying SBl 
orbits. The unpublished CfA orbits provided a sample of SBls that is approximately three 
times larger than that available in the literature, and that has a velocity reference frame 
which is not only self-consistent, but also consistent with our infrared frame. 



5. Discussion 

5.1. Validation of SB2 Solutions 

Four of the SB2 solutions shown in Figures [3] and H] depend on a single measurement of 
the secondary velocity, and another nine on only two measurements. Consequently, many 
of the derived and the resulting g's, have large uncertainties. Figure E] plots q^i^ 

against the measured g's to evaluate the plausibility of the measured values. As expected, 
q is greater than or within la of gmin for all of the binaries, except vB 59 and H382 which 
are consistent with ^minto better than 2a. Both vB 59 and H382 have long periods and 
were observed at orbital phases where the velocity separation was small, resulting in low 
precision measurements of K2. However, we have two observations of each system at slightly 
different phases and the measured secondary velocities are consistent. Therefore, we are 
confident that we are detecting the secondary in both of these systems. Observations of 
vB 59 in ~ 2012 and of H382 as soon as ~ 2009, when the orbital phases will have changed 
significantly, would further improve the SB2 solutions for these systems. 

The flux ratio, a, measured by our cross- correlation routine offers an additional valida- 
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tion of each SB2 solution. Figure [7] plots the a measured for each binary, averaged over all 
observations, against the measured mass ratio. Also shown are theoretical i/-band curves 
calculated from BCAH for binaries with primary masses from 0.6 Mq to 1.2 Mq. Two fac- 
tors complicate a direct comparison between our measured flux ratios and the theoretical 
values. First, the small wavelength range of CSHELL spectra severely limits our ability 
to measure precise flux ratios because individual spectral lines vary only a small amount 
relative to each other with changing spectral type. Considerably more accurate flux ratios 
can be measured with spectra covering a larger wavelength range, and thereby having many 
more spectral lines with differing dependencies on spectral type. Second, the theoretical 
curves represent the integrated flux over the entire i^-band, of which our CSHELL spectra 
only sample ~ 1.5%. Nonetheless, our measured values are well grouped along the curves. 
The two obvious outliers are L79 and vB 142. We have only a single, low S/N observa- 
tion of the long period binary L79. Our measured flux ratio is poorly constrained, varying 
significantly depending on the exact pair of templates used; the uncertainty shown in the 
figure probably underestimates the actual uncertainty. However, our measured secondary 
velocity is well constrained, independent of the template pair, and we are confident that the 
secondary detection is real. vB 142 has a very small mass ratio, making the measurement of 
the companion particularly difficult; §5.31 addresses this system in detail. 



5.2. Infrared Non-Detections 

Our infrared observations failed to detect the secondary in seven systems: vB 8, vB 30, 
vB 39, H411, L77, L90, and vB 114. Figure [8] shows the distribution of q-a^m for the infrared 
sample binaries; the systems not detected as SB2s are indicated by the hashed region. We 
expected that our sensitivity to binary companions would be incomplete for systems with 
the smallest mass ratios because these systems also have small flux ratios. However, several 
of the systems for which we did not detect the secondary have large gmin- For three addi- 
tional systems, vB 43, L57, and H509, we did not detect the secondary in a subset of our 
observations. We address each of these ten systems below. 

vB 8 and vB 30 have rapidly rotating F-type primary stars and both have small gmin, 
0.15 and 0.19, respectively. We observed both spectra in the infrared at multiple epochs, 
with orbital phases where the predicted velocity separation was large. We propose two 
possible explanations for our failure to detected these companions. The true value of q 
may actually be close to gmin- For a given mass ratio, the if-band flux ratio of a binary 
decreases as the mass of the primary increases. Because vB 8 and vB 30 have primaries 
more massive than those of a typical binary in our sample, their flux ratio's may be too 
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small to detect the secondaries at the S/N of our spectra. Alternatively, if the companions 
are also rotating rapidly, their spectral lines would be too broad to measure with CSHELL's 
limited wavelength coverage, even with high S/N observations. 

vB 114 and vB 39 have orbital periods of 4578 days and 5083 days, respectively. Our 
observations of both systems occurred at orbital phases such that the primary velocities were 
near 7, and the predicted velocity separations were only a few kms~^. Velocity measure- 
ments under such conditions are inherently difficult, and even had we detected the secondary 
components, the resulting K2^s would be poorly constrained. vB 114 will have a more favor- 
able orbital phase in ~2010; the phase of vB 39 will not improve until ~2013. Both have 
moderate gmin, 0.44 for vB 39 and 0.29 for vB 114, so the flux ratios should not impede in 
detecting these secondaries. 

H411, L77, and L90 fall towards the faint edge of the distribution shown in Fig. [H and 
our observations of these targets have S/N insufficient to detect their secondaries. H411 has a 
small Qmin, ~ 0.18, which corresponds to a minimum H-band flux ratio of only a few percent, 
and is at the limit of our best observation of this binary, with S/N~50. Because H411 is 
faint, obtaining a spectrum with CSHELL that has better S/N would require many hours of 
integration; such an observation could be carried out with a high resolution spectrometer at 
an 8 — 10 m telescope in a relatively small amount of time. L77 has a large gmin of ~ 0.61 
and we expected to be sensitive to its secondary. However, our single observation of L77 
has S/N~25, and occurred at an orbital phase when the velocity separation was small. The 
orbit of this system is such that even had we detected the secondary, the resulting SB2 orbit 
would be largely unconstrained. L77 will be at a more favorable orbital phase in 2009. L90 
has 

Qmin ~ 0.32; we observed it on two occasions, each with S/N~50. However, its long 
period and current orbital phase indicate a small velocity separation. This system will not 
be at a more favorable orbital phase until ~2011. 

We did not detect a companion in our observations of vB 43 on 2005 Sep 30, L57 on 
2005 Nov 25 and 2006 Feb 2, and H509 on 2005 Oct 2. Each of these observations were 
carried out at orbital phases corresponding to primary velocities near 7. The observation 
of vB 43 was compromised further by a low S/N, relative to our other observations of it. 
L57 has small Ki and K2, which when combined with it's K2 primary spectral type make 
it a difficult target to observe with CSHELL: future observations would benefit from higher 
spectral resolution or a larger free spectral range. 
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5.3. VB 142 

Our four infrared observations of vB 142 provided good phase coverage, and included 
measurements near the maximum velocity separation and on both sides of the 7-velocity. 
Despite this, the SB2 solution that we derive is poorly determined and the residuals, V2{fit) — 
V2{measured) , are large. We do not fully understand the reasons behind this poor solution, 
but we include it here for the following reasons. 

Our observations, except for that on 2005 November 28, return plausible velocities for 
the companion, albeit with large uncertainties. The 2005 November 28 observation occurred 
at phase ~ 0.40, and with small velocity separation, so our inability to accurately measure 
the secondary in this spectrum is not surprising. The average a that we measure for all of 
the observations, ~ 0.04, is, however, much larger than that predicted by BCAH (Figure [7]). 
To investigate this discrepancy we used our M-type template LHS2351 to introduce an 
additional spectrum into our observed vB 142 spectra from 2004 October 1 and 2005 February 
22. We added in this component with a "true" flux ratio, atme, ranging from 0.05 to 
0.0001, and, to avoid confusion with the actual vB 142 companion, at a radial velocity of 
— 20kms~^. We then attempted to recover this signal with our correlation procedure and 
the set of templates used in the original vB 142 analyses, excluding LHS2351. We recovered 
the LHS2351 spectrum at the proper velocity with atme as small as 0.0005. However, the 
uncertainty of the measured velocity increased as atme decreased, and the measured flux 
ratio, ameas. became unreliable for atme ^ 0.01. Consider, for example, the case of LHS2351 
introduced into the 2005 February 22 spectrum with a radial velocity of -20kms^^ and 
<^true = 0.005. Our correlation routine recovered this signal with a velocity of —25.1 ± 
7.0kms~^ and ameas. = 0.019 ±0.005. The larger than expected ameas. is consistent with the 
results obtained for the vB 142 companion, and may arise because our primary templates 
do not precisely match the vB 142 primary spectrum. For small flux ratios, the correlation 
routine tries to correct for this mismatch by scaling the primary using a; we have previously 



reported this behavior (IBender et al.ll2005l ). 



Our modeling with LHS2351 gives us confidence that we can detect a companion with 
a very small flux ratio. The q that we derive for vB 142, 0.06 ± 0.04, is consistent with 
9min ~ 0.04 from the SBl solution. Whatever the true value of q may be, our 3a upper limit 
of g < 0.18 is very small. Our primary goal in this endeavor is to measure the binary mass 
ratio distribution for the Hyades (§1): q for vB 142 is sufficiently well determined for this 
purpose. Of additional interest, the primary of vB 142 is a G5 star, and so a companion 
with q = 0.06 ± 0.04 could be a brown dwarf, which would be an important discovery in the 



Hyades (iGuenther et al.ll2005l ). 
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5.4. Component Masses 

While spectroscopic observations of an SB2 yield its dynamical mass ratio, they do not 
measure its orbital inclination and so alone they cannot provide a dynamical measurement 
of the individual component masses. Observations of the visual orbit measure the inclination 
and the total mass, and when combined with the mass ratio result in the individual masses. 
All of the binaries in our sample have components with a small angular separation and their 
visual orbits are not currently available. Those with periods longer than a few hundred days 
are resolvable with adaptive optics imaging at a large aperture telescope. 

In the absence of visual orbits, we can still obtain good estimates of the individual masses 
i f the distance is know n. Hipparcos measured the parallax of most of our sample binaries 



( iPerryman et al.lll998l ). and all have precise photometric measurements of their total flux 
from 2MASS at J, H, and K. We combined these with our measured mass ratios and a 
theoretical mass-luminosity isochrone from BCAH to calculate the individual component 
masses. We chose the BCAH models for several rea sons: they show a good, albeit not 



perfect, agreement with measured dynamical masses (IHillenbrandl l2004j ) : they include the 
effects of atmospheres, which are important in the low mass regime that applies to most of 
our secondaries; and lastly, they are provided in a convenient form that specifies magnitudes 
in the standard photometric bands used by observers. We used the 625 Myr isochrone, while 
noting that at such an old age the mass- luminosity relationship is mostly insensitive to age. 
Table [5] lists the calculated component masses. 

The uncertainties given in Table [5] include contributions from the parallax, photometry, 
and mass ratio; they do not include any uncertainties from the BCAH models. All of the 
2MASS J, H, and K photometric uncertainties are small, ~ 0.02 — 0.03 mag. Because our 
binaries have small flux ratios, the precision with which we determine the primary masses 
is mostly dependent on the precision of the parallax measurements. The uncertainty for the 
secondaries strongly depends on the precision of the mass ratios. Most of the primaries have 
masses determined to better than 10%, while some of the secondaries approach this level. 
Hipparcos did not measure the parallax of vB 59 or L57, so for these systems we used values 
reported in the Tycho catalog, and the resulting uncertainties on both the primary and 
secondary masses are large. The parallaxes of H69, H441, and L79 have not been measured, 
so we estimated their primary masses directly from their spectral type and assumed an 
uncertainty of ~ 0.1 Mq. The secondary masses then follow directly from our measured mass 
ratios. Finally, the primary of vB 68 is more massive than the ran ge covered by BCAH , so for 



this system only we used the empirical isochrone determined by iPinsonneault et al.l (120041 ) 



and note that the resulting masses are consistent with measured spectral types. 
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6. Summary 

We have obtained high resolution infrared spectroscopy of 32 SBs in the Hyades, whose 
SBl orbital parameters have been measured by RPS and DWL at the CfA, in order to 
detect their companions and thereby study the binary mass ratio distribution in this young 
cluster. We detected the companion in 25 of these systems. For these, we combined our 
results with the SBl parameters to determine their solutions as SB2s. Some of the SB2 
solutions we report have low precision for K2. However, obtaining precise orbital parameters 
for individual systems was not our objective here. Instead, our intent was to constrain the 
distribution of mass ratios in binaries with low mass companions, and our results are sufficient 
for this purpose. We also estimated the primary and secondary masses of our sample binaries 
using 2MASS photometry, Hipparcos parallax measurements, and our measured mass ratios. 
The mass ratios of the binaries with the most reliable SB2 solutions span the range from 
g ~ 0.1 — 0.8, corresponding to secondary masses as small as ~ 0.15 M©. We also detect a 
very low mass companion to vB 142. The solution for its mass ratio is not yet reliable, but 
it appears to be g < 0.18 at the 3a level, and may represent the detection of a brown dwarf 
companion. 

The precision of our derived primary masses is limited by the uncertainties in the Hippar- 
cos parallax measurements. The secondary mass measurements, however, can be improved 
significantly by reducing the uncertainties on the measured mass ratios through additional in- 
frared observations. Direct observations of sample members as visual binaries would measure 
their orbital inclinations and total masses, and when combined with the spectroscopy would 
yield dynamical component masses. By utilizing the Hipparcos distances, such measurements 



would contribute a test of the theoretical mass- luminosity relationships (e.g.. iMathieu et al. 



20071 ). Improving the SB2 solutions or measuring the visual orbits of our sample would 
require a significant commitment of observing time and analysis resources. The visual or- 
bit mapping may require technological improvements in interferometry and adaptive optics 
techniques. 

The results presented here are an initial step in a program to produce distributions of 
mass-ratio and secondary mass for Hyades cluster binaries with periods from a few days to 
a few thousand days. Future papers in this series will combine our new determinations of 
mass ratios with available orbital solution s from other Hyades b inary surveys, including the 



the spatially resolved systems studied by iPatience et al.l (Il998l ). to present the mass ratio 



and secondary mass distributions for the cluster. We also intend to address more fully the 
set of hierarchi c al trip le systems, for which, when combined with the speckle observations of 



Patience et al.l (119981 ). we have information on both the inner and outer orbits. 
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Table 1. Hyades Infrared Sample 



Target HD R.A. (J2000) Dec. (J2000) Sp. Type V H 



vB 8 


25102 


03 


59 


40. 


.49 


+10 


19 


49. 


.4 


F5 


6. 


.4 


5.4 


vB 9 




04 


00 


39, 


,54 


+20 


22 


49, 


.5 


G4 


8, 


.7 


7.0 


L20 


284163 


04 


11 


56, 


,22 


+23 


38 


10. 


.8 


KO 


9. 


.4 


6.6 


H69 




04 


12 


21, 


.44 


+16 


15 


03. 


.5 


Ml 


14 


.0 


9.1 


vB 30 


27397 


04 


19 


57, 


.70 


+14 


02 


06. 


.7 


FO 


5. 


.6 


4.9 


L33 


286770 


04 


22 


25, 


.69 


+11 


18 


20, 


.6 


K8 


9, 


.8 


7.1 


vB 40 


27691 


04 


22 


44, 


.17 


+15 


03 


21. 


.9 


GO 


7. 


.0 


5.6 


vB 39 


27685 


04 


22 


44, 


.78 


+16 


47 


27. 


.7 


G4 


7. 


.8 


6.3 


vB 43 


284414 


04 


23 


22, 


.85 


+19 


39 


31, 


,2 


K2 


9, 


,4 


7.3 


vB 59 


28034 


04 


26 


05, 


.86 


+15 


31 


27, 


.6 


G8 


7, 


.5 


6.2 


vB 62 


28033 


04 


26 


18, 


.50 


+21 


28 


13, 


.6 


F8 


7. 


,4 


6.1 


H382 


28068 


04 


26 


24, 


.61 


+16 


51 


12. 


.0 


Gl 


8. 


.0 


6.6 


H411 


285828 


04 


27 


25, 


.34 


+14 


15 


38, 


.5 


K2 


10 


.3 


7.8 


L57 


285766 


04 


27 


58, 


.96 


+18 


30 


00. 


.9 


K2 


10 


.2 


7.7 


vB 68 


28294 


04 


28 


23. 


.40 


+14 


44 


27. 


.5 


FO 


5. 


.9 


5.1 


vB 69 


28291 


04 


28 


37. 


.21 


+19 


44 


26. 


.5 


G5 


8. 


.6 


7.0 


H441 


285806 


04 


28 


50, 


.81 


+16 


17 


20, 


.3 


K7 


10 


.7 


7.6 


vB 77 


28394 


04 


29 


20. 


.55 


+17 


32 


41. 


.8 


F7 


7. 


.0 


5.8 


H509 


28634 


04 


31 


37. 


.10 


+17 


42 


35. 


.2 


K2 


9. 


.5 


7.3 


H532 


286839 


04 


32 


25, 


.65 


+13 


06 


47, 


.6 


KO 


11 


.0 


7.8 


vB 96 


285931 


04 


33 


58, 


.54 


+15 


09 


49, 


.0 


KO 


8, 


.5 


6.6 


L79 




04 


34 


10, 


.73 


+11 


33 


29. 


.6 


K7 


11 


.7 


8.3 


L77 




04 


34 


49, 


.76 


+20 


23 


41. 


.6 


K7 


11 


.1 


8.0 


vB 102 


29310 


04 


37 


31, 


.98 


+15 


08 


47, 


.2 


Gl 


7, 


.5 


6.1 


L90 


29896 


04 


43 


15, 


.70 


+17 


04 


08. 


.8 


KO 


9. 


.9 


7.5 


vB 142 


30246 


04 


46 


30. 


.39 


+15 


28 


19. 


.4 


G5 


8. 


.3 


6.8 


vB 113 


30311 


04 


46 


45. 


.58 


+09 


01 


02. 


.7 


F5 


7. 


.2 


5.9 


vB 114 


30355 


04 


47 


37, 


.57 


+18 


15 


31 


,4 


GO 


8, 


.5 


6.9 


vB 115 


284787 


04 


48 


42, 


.12 


+21 


06 


03. 


.6 


G5 


9. 


.1 


7.2 


vB 121 


30738 


04 


50 


48, 


.54 


+16 


12 


37. 


.6 


F8 


7. 


.3 


6.2 


vB 151 


240692 


05 


05 


40, 


.38 


+06 


27 


54. 


.6 


K2 


9 


.9 


7.6 


BD+02 1102 


40512 


05 


59 


29, 


.92 


+02 


28 


34. 


.2 


F5 


7. 


.8 


6.7 
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Table 2. Log of Infrared Observations 



Target 


UT Date 


JD-2400000 


Instrument 


V2 

( kms~^ 


) 


vB 8 


2004 Oct 4 


53282.967 


C 








2005 Sep 30 


53643.982 


C 






vB 9 


2004 Oct 2 


53280.978 


c 


31.8 ± 


2.1 




2004 Oct 4 


53282.993 


c 


32.9 ± 


1.8 


L20 


2004 Oct 1 


53279.959 


c 


133.9 ± 


4.0 




2004 Oct 2 


53280.950 


c 


- 53.6 ± 


3.6 




2004 Oct 3 


53282.148 


c 


115.6 ± 


4.3 




2005 Oct 10 


53645.985 


c 


84.6 ± 


9.7 




2005 Oct 3 


53647.112 


c 


-6.9± 


6.5 


H69 


2004 Oct 3 


53282.999 


c 


60.4 ± 


2.6 




2005 Feb 22 


53423.715 


N 


53.9 ± 


4.3 




2006 Feb 3 


53769.734 


c 


45.7 ± 


6.8 


vB 30 


2004 Oct 3 


53282.035 


c 








2004 Oct 4 


53283.036 


c 






L33 


2004 Oct 3 


53281.959 


c 


32.1 ± 


1.3 




2005 Oct 2 


53646.018 


c 


42.5 ± 


2.6 




2005 Nov 26 


53700.826 


c 


44.5 ± 


2.2 


vB 40 


2004 Oct 1 


53280.012 


c 


-45.5± 


1.9 




2004 Oct 3 


53282.123 


c 


132.9 ± 


3.0 




2005 Oct 2 


53645.958 


c 


117.7 ± 


5.8 




2005 Oct 3 


53647.020 


c 


52.9 ± 


4.6 


vB 39 


2005 Nov 26 


53700.897 


c 








2006 Feb 2 


53768.727 


c 






vB 43 


2004 Oct 3 


53282.052 


c 


42.2 ± 


2.5 




2005 Feb 22 


53423.754 


N 


48.0 ± 


2.0 




2005 Sep 30 


53644.008 


c 








2005 Nov 27 


53701.755 


c 


35.0 ± 


2.2 


vB 59 


2005 Nov 27 


53701.813 


c 


34.2 ± 


4.4 




2006 Feb 2 


53768.812 


c 


36.9 ± 


2.9 


vB 62 


2004 Oct 1 


53280.052 


c 


85.1 ± 


4.6 




2005 Nov 25 


53699.778 


c 


89.6 ± 10.9 




2005 Nov 28 


53702.758 


c 


-21.4± 


3.2 


H382 


2005 Nov 28 


53702.801 


c 


53.2 ± 


4.0 




2006 Feb 2 


53768.746 


c 


50.9 ± 


1.2 


H411 


2004 Oct 3 


53282.078 


c 








2005 Oct 2 


53646.060 


c 
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Table 2 — Continued 



Target 


UT Date 


JD-2400000 


Instrument 


V2 

( kms~^ 


) 


L57 


2004 


Oct 3 


53282, 


.104 


C 


46.7 ± 


2.0 




2005 


Nov 25 


53699 


.828 


c 








2006 


Feb 2 


53768, 


.889 


c 






vB 68*^ 


2005 


Oct 3 


53657, 


.074 


c 


22.8 ± 


1.1 


vB 69 


2004 


Oct 1 


53280, 


.082 


c 


-4.1± 


2.0 




2005 


Sep 30 


53644, 


.039 


c 


48.1 ± 


3.1 




2005 


Nov 27 


53701, 


.787 


c 


17.7 ± 


5.2 


H441 


2005 


Nov 28 


53702, 


.851 


c 


31.7± 


3.9 




2006 


Feb 2 


53768, 


.845 


c 


32.7 ± 


2.6 


vB 


2004 


Oct 1 


53280, 


.010 


c 


22.7 ± 


4.9 




2005 


Nov 30 


53644, 


.122 


c 


45.2 ± 


1.9 


H509 


2004 


Oct 2 


53281, 


.117 


c 


51.5 ± 


2.8 




2005 


Oct 2 


53646, 


.099 


c 








2006 


Feb 2 


53768, 


.786 


c 


27.0 ± 


2.0 


H532 


2004 


Oct 2 


53281, 


.015 


c 


-61.9± 


5.3 




2004 


Oct 4 


53283, 


.055 


c 


67.5 ± 


7.5 




2005 


Nov 26 


53700, 


.916 


c 


24.9 ± 10.2 


vB 96 


2006 


Feb 3 


53769, 


.922 


c 


44.6 ± 


1.6 


L79 


2006 


Feb 3 


53769, 


.792 


c 


45.9 ± 


2.8 


L77 


2006 


Feb 3 


53769, 


.840 


c 






vB 102 


2004 


Oct 2 


53281, 


.041 


c 


46.8 ± 


2.7 




2005 


Sep 30 


53644, 


.068 


c 


26.2 ± 


3.4 




2006 


Feb 2 


53768, 


.921 


c 


12.5 ± 


3.4 


L90 


2005 


Nov 25 


53699 


.881 


c 








2006 


Feb 3 


53769, 


.889 


c 






vB 142 


2004 


Oct 1 


53280, 


.118 


c 


54.2 ± 


2.1 




2005 


Feb 22 


53423, 


.860 


N 


17.2 ± 


6.4 




2005 


Oct 1 


53645, 


.099 


c 


29.0 ± 


3.5 




2005 


Nov 28 


53702, 


.910 


c 


20.2 ± 


1.8 


vB 113 


2004 


Oct 2 


53281, 


.066 


c 


50.7 ± 


1.4 




2005 


Oct 1 


53645, 


.061 


c 


51.5 ± 


3.4 


vB 114 


2005 


Nov 27 


53701, 


.930 


c 






vB 115 


2004 


Oct 4 


53283, 


.082 


c 


51.0 ± 


1.8 




2005 


Nov 27 


53701, 


.894 


c 


34.3 ± 


1.8 


vB 121 


2004 


Oct 1 


53280, 


.135 


c 


- 52.2 ± 


3.4 




2004 


Oct 3 


53282, 


.134 


c 


143.2 ± 


6.2 
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Table 2 — Continued 



V2 



Target 


UT Date 


JD-2400000 


Instrument 


(kms ^) 




2005 Nov 28 


53702.941 


C 


70.6 ± 12.7 


vB 151 


2004 Oct 4 


53283.112 


c 


28.7 ± 3.3 




2005 Oct 2 


53646.135 


c 


51.9 ± 1.9 




2005 Nov 25 


53699.929 


c 


54.9 ± 3.5 


BD+02 1102^ 


2004 Oct 4 


53283.136 


c 


24.5 ± 4.8 




2005 Oct 1 


53645.135 


c 


32.3 ± 6.2 



^Indicates a rapidly rotating primary. See §3.41 for a discussion. 



Table 3. Orbital Solutions 





r 








To 


7 


^1 




K2 








ai sin i 




02 sin i 


3 . 

iVfi sm % 




M2 


sin^ i 


larget 


(days) 


e 




;) (JD-2400000) (kms"! 


) \ Km s ) 


( km s 


) 




1 


(Gm 


) 




(Gm 


) 






(Mq) 


vB 9 


5070 


0.22 


107 




48747 


36.87 


2.74 


7.50 ± 


0.92 


0, 


,36 ± 0.05 


186 


± 


23 


510 




± 


84 


0.38 ±0.13 


0, 


,140 


± 0.034 


L20 


2.394358 


0.053 


286, 


,4 


50553.507 


40.86 


66.96 


98.2 


± 


3.6 


0, 


,68 ± 0.03 


2.202 


± 


0.014 


3, 


,23 ± 


0.12 


0.661 ± 0.053 


0, 


,451 


± 0.020 


H69 


128.114 


0.062 


337 




49110.9 


37.96 


13.80 


22.9 


± 


1.2 


0, 


,60 ± 0.03 


24.26 


± 


0.56 


40, 


,3 


± 


2.1 


0.407 ± 0.048 


0, 


,245 


± 0.019 


L33 


1044.9 


0.250 


295, 


,5 


48902 


40.11 


5.67 


8.16 ± 


0.58 


0, 


,69 ± 0.05 


78.9 


± 


2.6 


113, 


5 


± 


8.1 


0.153 ± 0.024 


0, 


106 


± 0.011 


vB 40 


4.000177 


0.0043 


118 




48353.24 


37.920 


39.88 


87.5 


± 


2.7 


0, 


,46 ± 0.02 


2.1934 


± 


0.0042 


4, 


,81 ± 


0.15 


0.588 ± 0.043 


0, 


,268 


±0.011 


vB 43 


589.76 


0.619 


304. 


,7 


50002.65 


38.925 


9.61 


14.5 


± 


1.2 


0, 


,66 ± 0.05 


61.23 


± 


0.78 


92, 


,4 


± 


7.4 


0.249 ± 0.046 


0. 


,165 


±0.017 


vB 59 


5724 


0.975 


224. 


,2 


51385 


39.349 


15.1 


84 




28 


0, 


.18 ± 0.06 


264 


± 


28 


1460 




± 


490 


5.4 ± 4.8 


0. 


,96 


±0.56 


vB 62 


8.550647 


0.212 


41, 


,0 


49925.248 


38.247 


16.73 


72.6 


± 


8.8 


0, 


,23 ± 0.03 


1.923 


± 


0.017 


8, 


,3 


± 


1.0 


0.48 ±0.15 


0, 


,110 


± 0.022 


H382 


2657 


0.682 


248, 


,4 


48884.7 


38.73 


7.79 


21.1 


± 


2.4 


0, 


,37 ±0.05 


208 


± 


11 


563 




± 


66 


1.90 ±0.55 


0, 


,70 


±0.14 


L57 


1911.3 


0.486 


94, 


,1 


49211.6 


39.124 


6.63 


7.7 


± 


2.1 


0, 


,86 ± 0.23 


152.2 


± 


2.6 


177 




± 


48 


0.21 ±0.21 


0, 


180 


± 0.053 


vB 68 


331.66 


0.288 


319 




50290.4 


40.00 


11.19 


16.1 


± 


2.4 


0, 


,70 ± 0.11 


48.9 


± 


2.5 


70 




± 


11 


0.36 ±0.12 


0, 


,251 


± 0.051 


vB 69 


41.6729 


0.643 


328, 


,4 


49443.176 


38.993 


7.03 


44.0 


± 


5.9 


0, 


,16 ± 0.02 


3.08 


± 


0.11 


19, 


,3 


± 


2.6 


0.222 ± 0.082 


0, 


,0355 


± 0.0085 


H441 


7494 


0.186 


324 




52730 


40.44 


3.62 


8.0 


± 


0.3 


0, 


,45 ± 0.03 


366 


± 


1 T 

17 


807 




± 


36 


0.796 ± 0.073 


0, 


,360 


± 0.033 


vB 77 


238.86 


0.200 


132 




48557.3 


39.22 


6.58 


20.0 


± 


1.8 


0, 


,33 ±0.03 


21.18 


± 


0.91 


64, 


,3 


± 


5.7 


0.329 ±0.075 


0, 


108 


± 0.016 


H509 


849.95 


0.174 


316, 


,9 


48643 


39.553 


6.35 


13.7 


± 


4.0 


0, 


,46 ±0.14 


73.0 


± 


1.5 


157 




± 


46 


0.46 ±0.32 


0, 


,215 


± 0.086 t 


H532 


1.484698 


0.0035 


267 




48170.37 


40.31 


68.94 


94 


± 


15 


0, 


,73 ±0.12 


1.4075 


± 


0.0044 


1, 


,92 


± 


0.30 


0.38 ±0.13 


0, 


,281 


± 0.052 


vB 96 


5100 


0.664 


310, 


,4 


50106 


41.14 


4.77 


9.6 


± 


3.5 


0, 


,50 ±0.18 


249.8 


± 


8.1 


500 




± 


180 


0.44 ±0.37 


0, 


,22 


±0.11 


L79 


3688 


0.765 


341. 


,8 


48662.1 


41.454 


4.31 


27 


± 


19 


0, 


,16 ±0.11 


140.8 


± 


7.9 


890 




± 


620 


2.7 ±5.2 


0. 


,43 


±0.53 


vB 102 


734.79 


0.513 


337. 


,2 


50015.8 


40.192 


3.92 


31.1 


± 


5.9 


0, 


,13 ±0.02 


34.0 


± 


1.4 


270 




± 


51 


1.84 ±0.97 


0. 


,231 


± 0.080 


vB 142 


975.7 


0.675 


271, 


,9 


50383.6 


41.684 


1.19 


21 


± 


14 


0, 


,06 ± 0.04 


11.8 


± 


1.3 


220 




± 


140 


0.42 ±0.82 


0, 


,024 


± 0.031 


vB 113 


2429 


0.327 


294, 


1 


51404 


41.544 


3.36 


11.1 


± 


2.9 


0, 


,30 ±0.08 


106.1 


± 


4.3 


350 




± 


91 


0.49 ±0.32 


0, 


150 


± 0.060 


vB 115 


1208.2 


0.480 


129, 


1 


50731.6 


40.730 


5.65 


10.9 


± 


3.7 


0, 


,52 ±0.18 


82.3 


± 


1.5 


158 




± 


54 


0.25 ±0.20 


0, 


131 


± 0.059 


vB 121 


5.750872 


0.361 


42, 


,0 


49743.257 


41.28 


20.11 


122.4 


± 


9.2 


0, 


,16 ±0.02 


1.483 


± 


0.015 


9, 


,03 ± 


0.68 


1.20 ±0.25 


0, 


197 


± 0.026 


vB 151 


629.37 


0.297 


203, 


,8 


50544.0 


42.167 


5.22 


9.8 


± 


2.2 


0, 


,53 ±0.12 


43.14 


± 


0.89 


81 




± 


18 


0.125 ±0.065 


0, 


,067 


± 0.020 


BD+02 1102 


32.5121 


0.344 


228, 


,7 


51056.33 


47.60 


21.31 


53 


± 


21 


0, 


,40 ±0.15 


8.94 


± 


0.21 


22, 


,4 


± 


8.6 


0.82 ±0.78 


0. 


,33 


±0.19 
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Table 4. Alternate Orbital Solutions 





P 


LO 


To 


7 




K2 




Target 


(days) 


e (deg) 


(JD-2400000) (kms-i) 


( kms"-*^) 


(kms^"'^) q 


SBl Rcf. 


L20 


2.394357 0.057 279.2 


43892.36 


36.76 


66.2 


98.8 ± 2.7 0.67 ±0.02 


G81 


vB 40 


4.00050 


0.060 13 


22274.81 


37.4 


36.1 


88.2 ± 2.7 0.41 ±0.02 


S21 


vB 43 


590.6 


0.638 303.1 


43512.9 


39.81 


9.91 


14.1 ± 1.3 0.70 ±0.07 


G85 


vB 62 


8.55089 


0.233 38.0 


42588.22 


38.77 


16.46 


87 ±30 0.19 ±0.06 


G78 


L57 


1907 


0.485 95 


43470.5 


40.23 


6.83 


6.8 ± 2.2 1.00 ±0.32 


G85 


vB 69 


41.6625 


0.662 326.9 


43650.67 


39.81 


7.28 


69 ±17 0.10 ±0.03 


G85 


vB 77 


238.87 


0.242 127 


43298 


39.81 


6.53 


20.1 ± 1.9 0.32 ±0.03 


G85 


H509 


844.6 


0.148 325 


44413 


40.32 


6.20 


12.6 ± 1.6 0.49 ±0.06 


G85 


vB 121 


5.75096 


0.354 54.9 


42192.06 


42.74 


19.70 


126.8 ± 8.1 0.16 ±0.01 


G78 



References. — G78 -ICriffin fc Gun^ |l978l ): G81 - iGriffin fc GuniJ |l98ll ): G85 - iGriffin et al 
|l985|); S21 - ISanfordI (|l92ll ) 
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Table 5. SB2 Component Masses 



Target Mi ( Mq ) M2 ( Mq ) 



vB 9 


1.15±0 


05 


0.41 ±0 


06 


L20 


0.88 ±0 


04 


0.60 ±0 


04 


H69^ 


0.45 ±0 


10 


0.27 ±0 


06 


L33 


0.93 ±0 


05 


0.64 ±0 


06 


vB 40 


1.38 ±0 


10 


0.63 ±0 


05 


vB 43 


1.01 ±0 


06 


0.67 ±0 


06 


vB 59^ 


0.95 ±0 


08 


0.17 ±0 


06 


vB 62 


1.23 ±0 


03 


0.28 ±0 


04 


H382 


1.09 ±0 


04 


0.40 ±0 


06 


L57'^ 


0.77 ±0 


09 


0.66 ±0 


19 


vB 68^^ 


1.70 ±0 


08 


1.19±0 


20 


vB 69 


1.03 ±0 


04 


0.16 ±0 


02 


H441'^ 


0.61 ±0 


10 


0.27 ±0 


05 


vB 77 


1.29 ±0 


03 


0.43 ±0 


04 


H509 


0.83 ±0 


04 


0.38 ±0 


12 


H532 


0.79 ±0 


08 


0.58 ±0 


15 


vB 96 


1.00 ±0 


05 


0.49 ±0 


18 


L79^ 


0.61 ±0 


10 


O.lOiO 


07 


vB 102 


1.19±0 


04 


0.15±0 


04 


vB 142 


l.lOiO 


04 


0.07 ±0 


04 


vB 113 


1.17±0 


03 


0.32 ±0 


12 


vB 115 


0.97 ±0 


05 


0.50 ±0 


18 


vB 121 


1.31 ±0 


04 


0.21 ±0 


03 


vB 151 


1.03 ±0 


09 


0.55 ±0 


13 


BD+02 1102 


1.29 ±0 


05 


0.52 ±0 


20 



^Estimated from primary spectral type 
"^Tycho parallax 

HPinsonneault et al.l ([200J) isochrone 
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Fig. 1. — The distribution of 2MASS //-magnitudes for the infrared sample binaries. The 
hashed region indicates binaries for which we did not detected a companion ( §5.21) . 
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^1 (predicted-measured) S ) 



Fig. 2. — The distribution of the difference between the primary velocity predicted from the 
visible SBl solutions and that measured from the infrared spectra. The outlying measure- 
ments with large velocity difference are attributable to spectra with low S/N. The distribution 
has a mean value of ~ 0.3 km s^^, and is fit by a Gaussian with with ~ 0.9 km s^^, indicating 
that the two reference frames agree to better than the infrared measurement uncertainties. 
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Fig. 3. — Double- lined velocity versus phase curves for the Hyades infrared sample, from R.A. 

4'»3i"^. The circles indicate the secondary velocities and associated uncertainties 
measured in the infrared. For each binary, the dashed curve shows the CfA SBl solution; 
the solid curve shows the SB2 solution. 
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Fig. 4.— Same as Figure [3l from R.A. 4'*32" to 5''59'" 
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Fig. 5. — (a) The distribution of mass ratios measured in the infrared for the SB2s in 
TableO (b) Same as (a), but includes the large variation in the uncertainties of our measured 
mass ratios. We distributed each mass ratio over a Gaussian with width equal to the la 
measurement uncertainty, clipped at q = gmin — 2 x cTq^^.^^ and q = I, and normalized to unity. 
The uncertainties shown for both (a) and (6) are 
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Fig. 6. — The minimum mass ratio, determined from f{M), plotted against the measured 
mass ratio. The dotted hne indicates gmin = q. 
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Fig. 7. — The measured flux ratio plotte d against the measur ed mass ratio. The curves show 
the theoretical if-band flux ratios from iBaraffe et al.l (1l998l ) for 625 Myr old binaries with 
primary masses of 0.6 Mq [solid), 0.8 Mq (dotted), 1.0 Mq [dashed), and 1.2 MQ{dot- dashed) . 
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Fig. 8. — The distribution of gmin for the systems detected (open) and not detected (hashed) 
as SB2s in the infrared. 



